A two-atom electron pump 
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The fabrication of single atom transistors [1- 
[4] paved the way for electronics based on single 
dopants [5, 6j. Recently the spectrum of a single 
dopant was measured electrically by coupling two 
such devices [7] . The next step towards promising 
functionalities for future nanoelectronics [8j con- 
sists in manipulating a single electron over two 
dopants. Here we demonstrate electron pump- 
ing through two phosphorus donors in series im- 
planted in a silicon nanowire. While quantized 
pumping is achieved in the low frequency adia- 
batic regime, we observe remarkable features at 
higher frequency when the charge transfer is lim- 
ited by the different tunneling rates. The tran- 
sitions between quantum states are modeled in- 
volving a Landau-Zener transition, allowing to re- 
produce in detail the characteristic signatures ob- 
served in the non-adiabatic regime. 

Adiabatic charge pumping at zero d.c. bias, in contrast 
to stationary electronic transport at finite bias, results in 
a net direct current I = ef when exactly one electron is 
transferred during a cycle of period j (Refs. |9HTT|). Mo- 
tivated by metrological applications [12], there has been 
a growing interest in the non-adiabatic regime where this 
period is comparable to the characteristic time which par- 
ticles spend in the system p!3HT5] . All experiments so far 
were in either one of these regimes but were not subse- 
quently tuned to both. The ultra-scaled electron pump 
presented here, where the usual large dots are replaced 
by real nanometer-size atomic orbitals, offers this unique 
opportunity 

Single atom transistors are now realized in differ- 
ent ways, including both random [16 or determinis- 
tic m (TTJ placement of the dopant atoms. Electron 
pumping through many donors has been reported in a 
larger structure [18]. Here, despite the stochastic implan- 
tation process, transport through two dopants is obtained 
because 3 gates control an extremely small volume of sil- 
icon (60x20x60 nm^) in which the number of P dopants 
is only a few tens for a concentration of lO^^cm"^, be- 
low the metal-insulator transition (see Fig.jl]). The parts 
of the wire uncovered by the gates and adjacent spac- 
ers are doped with arsenic donors at high concentration 
(^ 10^^ cm~^) to define conducting source and drain 
electrodes, separated by 60 nm. This number prevents 
transport through a single donor, which requires a much 



shorter structure [2 , but favors transport through two 
donors in series. Indeed a distance of ~ 30 nm between 
the two donors yields a mean tunnel coupling of the or- 
der of 1 GHz [19 , a good compromise between weak cou- 
pling and measurable current. Moreover we obtain sepa- 
rate control over the energy levels of the two dopants by 
splitting the top gate into two, which face each other at 
a distance of only 15 nm. Finally a voltage V^g is applied 
to the substrate [7 in order to tune the coupling between 
the dopants and the leads and rule out rectification ef- 
fects [20] , measured for instance in Ref. [211 Starting from 
a nearly perfect yield in terms of transistor behaviour, the 
final yield for two-dopant features remains higher than 
50% for the geometry and doping level chosen here. 

The sample is measured at 150 mK in a dilution re- 
frigerator. The source-drain current / is first measured 
with a fixed source-drain bias voltage Vd to select two 
ideally placed dopants. The back gate voltage is changed 
until a set of triangles, characteristic of two levels in se- 
ries [22 , is obtained when scanning the front gates (see 
Fig. [l]i). This coupled donor system allows to perform 
atomic spectroscopy, in contrast with a single donor tran- 
sistor based on the same technology where fluctuations 
of the local density of states dramatically affect the cur- 
rent at finite bias [23|. A large (lOmeV) separation in 
energy between the ground and first excited states was 
measured in the same sample [7] as investigated here. 
Besides this large valley-orbit splitting which makes it 
quite insensitive to a small bias, our 2-donor system is 
also well isolated in energy from any other states which 
could carry a current from source to drain (see Fig. [l]i). 
The two triple points are separated by a Coulomb repul- 
sion energy Eq = 4.0 ± 0.1 meV, which is markedly larger 
than any value reported so far [2T. 

We now set Vd to zero and modulate the top gates 
(i = 1, 2) with phase shifted sinewaves added to the d.c. 
voltages Fg, = V^^ + V^^ sin(27r/t ± (p/2), in order to 
follow elliptical trajectories in the {Vgi, Vg2) plane en- 
closing a triple point [^(see Fig. [l^). 

Adiabatic pumping occurs when the timescale of the 
driving is much larger than the timescale associated to 
tunneling to the leads for an electron in the 2-donor sys- 
tem, <C 1, and when there are no excitations in 
the internal dynamics due to the driving signal, requir- 
ing 2iTt\2/hv ^ 1 (see Eq. [l]). Here F is a tunneling 
rate, 5e = ei — €2 the difference in energy between the 
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FIG. 1: Control of two dopants with three gates, a, Schematic tilted view of the sample along the source-drain axis, 
with 2 phosphorus dopants in yellow in the small silicon region controlled by 3 gates, b, Transmission electron micrograph of 
the cross section of a nominally identical sample. The two top gates and the wire are highlighted by blue and white dashed 
lines respectively. The cross section is slightly offset from the gate region (along the red line in a) where the wire is thicker 
than below the gates because of the source/drain raised by epitaxy. The wafer substrate, separated by a 145 nm thick Si02 
layer, is used as a back gate, c. Schematic top view of the sample including the gate and source-drain control parameters, 
d. Source-drain current (/) with Vd = 2mV, Vbg = 11.5 V and superimposed axes giving the energy of each ionized donor with 
respect to the Fermi level in the leads. The two triangular triple points of the two-dopant system are isolated from any other 
visible single electron feature by at least 20 meV. e. Triple points at small bias (Vd = 50 /iV) measured with a lock-in amplifier. 
The black dashed lines indicate regions of stable charge states (ni, 712), where rii is the number of electron on the i*^ donor. 
The white dashed line shows an example of pumping trajectory realized in Fig. [S]:, around the working point indicated by the 
black cross. 



two levels, v = \ ^\ the speed at which they cross, ksT 
the energy scale set by the temperature and ti2 the tun- 
nel coupling between the two levels. Then the system 
follows the charge states given by the electrostatic sta- 
bility diagram (Fig.[l^). Fig. [2] shows the normalized d.c. 
current I/ef measured in this regime. As expected we 
observe two circular regions of constant quantized cur- 
rent / = ±e/ wherever a triple point is enclosed by the 
ellipse, and / = elsewhere. The sign of the current is 
reversed by changing the direction of the trajectory. The 
2- atom system offers intriguing advantages compared to 
other quantized electron pumps, due to the large energies 
at play in this system. The large Coulomb repulsion en- 
ergy between the dopants ensures wide and flat plateaus 
for the pumped current versus Vd'. a slope of 1.9 GQ is 
obtained over a Vd range of 3.1 mV at 1 MHz (see Supple- 



mentary Information, Fig. SI). In contrast with metal- 
lic Coulomb islands co-tunneling events via intermediate 
virtual states are strongly suppressed due to the 10 meV 
single level-spacing for each of the dopants [7 . Both these 
elastic cotunneling effects as well as remaining unwanted 
events due to inelastic tunneling are shown to be small. 

At higher frequency the pumped current deviates from 
the fundamental relation / = e/. Figs. [3^,c show the cur- 
rent map for / = 10MHz ^0\. The two circular regions 
appearing when enclosing the triple points are still visi- 
ble in Fig. |3]3, but the current inside is not independent 
of the working point anymore. This is the first manifes- 
tation of the deviation from the adiabatic regime. The 
second remarkable feature is the occurence of extra cur- 
rent outside these circular regions. This effect becomes 
dramatic in Fig. [3^, where, due to the vanishingly small 
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FIG. 2: Adiabatic electron pumping though the two 
dopants, a, Map of the pumped current at 1 MHz and zero 
bias voltage Vd. b, Pumped current for different frequencies 
along the Vg path indicated by the arrow in a. Quantized 
pumping / = zbe/ is indicated by horizontal dashed lines. 



surface area of the pumping cycle for cp = 175°, no signal 
is expected in the adiabatic case. 

The complex features observed in the non-adiabatic 
regime can be fully understood by means of the model 
depicted in Fig. [3^ and described in detail in the Supple- 
mentary Information. The transitions occuring within 
the two dopants are described through the dynamics of 
the closed two-level system. A dissipative part of the dy- 
namics accounts for the tunnel coupling to the source and 
drain. The transition rates between donor states due to 
tunneling to the leads are obtained from Fermi's golden 
rule [25 and are proportional to the tunneling rates to 
source and drain (Fs and Fd). Hopping of an electron 
between the two dopants, with tunnel amplitude ti2, can 
occur when the two levels cross each other and the result- 
ing time evolution of the two-level system is described ac- 
cording to the Landau-Zener problem [26 . Consequently, 
the probability of staying on the same dopant after the 
crossing, i.e. the transition from ground to excited state, 
is given by the Landau-Zener transition probability Plz- 
Then the probability for an electron to move from one 
dopant to the other when the levels are aligned depends 
on the speed v and is given by, 

l-PLZ = l-exp(^-2^|^^ (1) 

The strict separation of the dynamics of the closed system 
and of the coupling to the reservoirs is possible because 
the time during which the levels cross, ~ 2ti2/v^ is of the 
order of 0.1 ns and thus much smaller than the time scale 
related to the dissipative dynamics, Fg^,F^^. 

Using the electrostatic couplings obtained from the 
d.c. measurements we reproduce the various features ob- 
served in the experiment with a unique set of parame- 



ters (Fd = 670MHz, Fs = 18MHz and tis/Zi = 610 MHz), 
as illustrated in Figs. [3]3,d for (^ = 175° and (p = 100°. 
The whole set of phase differences (p is simulated with 
the same fidelity (see Supplementary Material, Anima- 
tion 1). If nonadiabatic, the system does not necessarily 
relax to the ground state, either due to slow tunneling to 
the leads or due to its internal dynamics (Plz differs sig- 
nificantly from zero, if the condition 27rt^2/^'^ ^ 1 is no 
longer fulfilled). Each case yields characteristic features 
we discuss in the following. 

We focus on the working point indicated by the black 
cross in Fig. [3]i. In this region no current is detected 
in the adiabatic regime. The time evolution of the en- 
ergy levels is shown in Fig. [sjF fSV. Dopant 2 can ac- 
cept electrons only from the drain and not from dopant 
1 because it always lies below the Fermi level. After the 
electron moves to the dopant 1, its energy can be raised 
above the Fermi level, and it is eventually released into 
the drain lead. The key point here is that the transi- 
tion probability between the two discrete dopant levels 
is different from one, otherwise the electron sitting on 
dopant 1 would go back to the initial dopant 2 after the 
second level crossing. This missed tunnel event due to 
the Landau-Zener transition occuring with probability 
Plz 7^ in the non-adiabatic regime increases signifi- 
cantly the error rate [27j. Indeed, for this working point, 
we get Plz ~0.25 from the simulation for both cross- 
ings. In contrast, at low enough frequencies the electron 
always stays on the lowest energy level according to the 
adiabatic theorem (Plz ^ 0)- Also the rightmost red 
non-adiabatic feature can be understood in these terms. 

The non adiabatic feature inside the circle is of dif- 
ferent origin: here, it is the tunneling dynamics which 
matters. Even though 62 is driven from above to below 
the Fermi energy before ei, the slow tunneling dynamics, 
27rfv ~ F^/cbP, allows to charge dopant 1 first instead of 
2. This effect is important in particular if the coupling 
between dopant 1 and the drain is much stronger than 
the coupling between dopant 2 and the source. Indeed, 
Fd ^ Fs is fulfilled in the experiment, giving rise to the 
sign change within the circle. Note that this effect would 
survive even in the limit of Plz ~ 0, hence it does not 
rely on a non-adiabatic level crossing. The additional fea- 
tures can be understood analogously and are influenced 
by the individual coupling strengths Fs and Fd, with 
the exception of one feature that will be discussed in the 
following. 

In order to reproduce the faint side triangle along the 
anti-diagonal, visible in Fig. |3]i, we finally add a small 
effective inelastic relaxation rate Fin = 14 MHz between 
levels to the model. However this rate has a very limited 
impact on the overall pumped current. For the working 
point indicated by the cross, the simulation shows that 
only 10 % of the total current comes from this inelastic 
term, the remaining 90 % being due to the elastic cou- 
pling ti2' A much higher inelastic rate would suppress 
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FIG. 3: Non-adiabatic electron pumping at 10 MHz: experiment and simulation. a,c, 2D map of the current 
measured with — 175° (a) and 100° (c) between the two driving signals as a function of the working point (Vg^^, V^^). 
b,d, Simulation of the current with Fd = 670 MHz, Fs = 18 MHz and ti^jh — 610 MHz. The shape of the trajectory in the (V^i, 
V^2) plane is shown at the top left of the simulation for each phase difference (not to scale), and triple points are indicated with 
a dot. e. Sketch of the model used for the previous simulation. Tunneling to the leads occurs with rates Fd and Fs. Coupling 
between the two dopants by a tunnel amplitude ti2, leads to a Landau-Zener transition when the two levels cross each other, f. 
Calculated energy levels as a function of time for the working point indicated by the black cross in d. Landau-Zener transitions 
occur as immediate events when the two levels cross (between the dashed lines). 



the current coming from the Landau-Zener transition. 

Importantly, the ultimate atomic size of the sys- 
tem is crucial to get large level separation and strong 
Coulomb interaction. This enables us to use driving am- 
plitudes corresponding to an energy level modulation of 

2.5 me V, much higher than the electronic tempera- 
ture (300 mK) in the contact, equivalent to 25 /ieV. As 
a consequence the pumped current signal has remark- 
ably sharp borders between different pumping regions. 



allowing a detailed study of the mechanisms at play in 
each of these regions. Moreover, due to the large ampli- 
tudes, inter-level transitions occur far away (with respect 
to /cbT) from the resonances with the leads, which jus- 
tifies the zero temperature approach for our model. Fi- 
nally, since the crucial parameter determining the non- 
adiabaticity is the speed of the crossing, the actual fre- 
quency to observe these effects depends on the size of the 
ellipse. With the large amplitudes we use the crossover 
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is significantly lowered towards lower frequencies, making 
it easily accessible experimentally. 

We have investigated electron pumping through two 
dopant atoms connected in series and addressed with 
three gates. We achieved quantized pumping in the adi- 
abatic regime and, studying the crossover to the non- 
adiabatic regime, observed additional features in the 
pumping current. A model taking into account the tun- 
nel couplings to the leads and between the two donors 
energy levels reproduces all the details of the measured 
pumped current. This ultimate size system behaves like 
a nearly perfect two-level system, and the full control of 
the well separated states in the dopants opens new op- 
portunities for the study of pure spin pumping in a quan- 
tum system contacted to ferromagnetic leads [28] and for 
transporting quantum information [29] . 
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